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Abstract. Analytic energy bounds for Af-boson systems governed by 
semirelativistic Hamiltonians of the form H = X/i— i^Pi ~'~ ™^)^''^ ~ 
^2^—i<j^/^^j' V > 0, are derived by use of Jacobi relative coordinates. 

For gravity ii = c/N, these bounds are substantially tighter than earlier bounds 
and they are shown to coincide with known results in the nonrelativistic limit. 
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1. Introduction: the A^-body problem 

One-body Hamiltonians H composed of the relativistic expression \/p'^ + for the 
kinetic energy of particles of mass m and momentum p and of a coordinate-dependent 
static interaction potential V'(r), defined as operator sum 

H = v/p2 + m2 + y(r), 

provide a simple but very efficient tool for the description of relativistically moving 
particles 12 El E El E) • They have been used, for instance, for the description 
of hadrons as bound states of quarks [3 El- One of the advantages of this kind 
of semirelativistic treatment is that its generalization to the many-body problem 
is straightforward to formulate. A semirelativistic Hamiltonian for a system of N 
identical particles interacting by pair potentials V{rij) is given by 

N N 

H = Y.^pj+m^+ V{u,). (1) 

i—l l—i<j 

We shall use the notational simplification p = |p|, r = |r|, orr^j = |ri—rj|, whenever no 
ambiguity is introduced by doing so. Many approaches to such many-body problems 
for identical particles employ the very powerful constraint of permutation symmetry to 
generate their reduction to a two-body problem with a Hamiltonian Ti. whose spectrum 
is used to approximate the many-body energy eigenvalues or to generate a lower 
energy bound. This reduction may be effected in various ways, which leads to the 
problem of finding the most effective reduced problem, the one which would provide 
the highest lower bound. In this paper, we use Jacobi relative coordinates which 
introduce technical difficulties but yield a lower bound that is much improved over 
previous results. We have already achieved a similar improvement for the harmonic- 
oscillator potential 9 and for convex transformations of it ^01. Here, we show that 
one can derive tight A^-body energy bounds for the (more physically important) case of 
the gravitational potential: V{r) = —v/r, v > 0. There is a long history of attention 
to the corresponding nonrelativistic problem, dating back at least to 1967 I12|. 
Also the semirelativistic case has been discussed before, notably by Lieb and Thirring 
|13|. by Lieb and Yau (TJI, as well as by Martin and Roy 15 . 

In their rigorous investigation JJl of Chandrasekhar's theory of the stellar collapse 
and its formal connection to quantum physics, Lieb and Yau were able to demonstrate 
exactly that in the simultaneous limit N oo and u — > of particle number N and 
gravitational coupling constant v such that their product Nv is kept fixed (at an 
arbitrary value below some critical value) the Schrodinger equation reduces to the 
"semiclassical" approximation represented by a Hartree-type equation for the density 
of the bosons. In particular, they succeeded in proving the convergence of the lowest 
quantum energy, E, defined as the infimum of the spectrum of the Hamiltonian H given 
by Eq. to the corresponding semiclassical minimum Hartree energy. Interestingly, 
precisely such kind of behavior in the limit N oo has been conjectured already 
before in Ref. 

The study in Ref. JSl shows that for the one-body case with a gravitational 
potential a lower bound to the energy spectrum may be expressed in terms of 
the lowest energy of a Schrodinger operator with a Kratzer potential, which is a 
Coulomb potential "spiked" with an additional term of the form A/r'^. In spite of the 
complications of relative coordinates but with boson permutation symmetry in the 
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individual-particle coordinates, we may take advantage of this result to construct 
improved lower energy bounds for systems comprised of identical particles with 
integer spin. 

The plan of this paper is as follows. After recalling, in Sec. 13 existing Coulomb 
one-body lower energy bounds, and analyzing, in Sec. |21 the simple {N/2) N- 
body lower energy bound, we reformulate, in Sec. 01 the A^-particle problem in 
relative coordinates. This allows us to derive, in Sec. an improved lower energy 
bound. To this result we adjoin, in Sec. El a variational (Rayleigh-Ritz) upper energy 
bound obtained with the help of a scale-optimized Gaussian trial wave function. We 
summarize our results in Sec. |3 and inspect, in particular, their large- TV limit. 



2. One-body lower energy bounds for the relativistic Coulomb problem 

The case of a one-particle Hamiltonian with a Coulomb potential V{r) = —v/r has 
been extensively investigated. We know from the pioneering work of Herbst that 
the Hamiltonian is self-adjoint for v < and has a Friedrichs extension up to the 
critical value v — 2/tt. Herbst obtains the following lower bound to the ground-state 
energy E: 

For the smaller range < < ^ of the Coulomb coupling v, this result was 
strengthened by Martin and Roy ^H] to 



E>., f^^^ , 0<.<1, (2, 

This energy bound was found by considering the squares of the kinetic-energy and the 
potential-energy terms and by relating certain expectations to the previously studied 
exact solution of the Klein-Gordon Schrodinger equation. Let us now briefly re-derive 
this same result, emphasizing that this lower bound is, in fact, generated by the lowest 
energy of a Schrodinger operator JJ]; this latter problem had earlier been analyzed 
by Kratzer and others [H EOl Ell 122 ESI El . 

We suppose that the exact normalized ground state of is "0 so that {p^ + 
jji2^i/2^ = (i? + v/r)tp. By considering the equality between the squares of these 
equal vectors we obtain, upon introducing a quantity £ by £{£ + 1) = — w^. 



the inequality on the right-hand side of Q arises from the variational principle applied 
to the Kratzer Hamiltonian: 

The well-known expression for the bottom, for given of the hydrogen energy 
spectrum remains valid for negative non- integer £ provided \£\ < i, which is equivalent 
to the constraint w < i. By solving for £ in terms of the coupling parameter v, we 
recover the Martin-Roy lower bound (EJ. 
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3. "Simple" (or N/2) lower energy bound for general A^-body problems 



This rather simple lower energy bound is not limited to the gravitational pair potential, 
and it allows us to prove that a given A^-body Hamiltonian is bounded from below. 
Applying the same reasoning to the (soluble) Schrodinger harmonic-oscillator problem, 
defined by the Hamiltonian 



N 



H 



N 

i=l l—i<j 



E 



with the exact ground-state energy E — 3{N — 1)^/N, one gets only E/ \/2 whereas a 
general lower bound based on Jacobi relative coordinates yields indeed the exact energy 
for this potential ■ A lower bound to the ground-state energy of the Hamiltonian 
is provided by the bottom £n/2 of the spectrum of the one-body Hamiltonian 
operator 

— ^ A^- 1 



n 



N/2 



N 



pi 



-V{r) 



(4) 



since boson permutation symmetry of the exact TV-body normalized wave function 5* 
implies (^', H'^) = (5*, h'^), where the operator h is a. two-body Hamiltonian given by 



pf + m2 + Jpl + {N - l)V{ri2) 



Changing the coordinates of this two-particle problem to r = ri — r2 and R = ri + r2, 
the individual momenta are given, in terms of the corresponding total and relative 
momenta variables, by pi^2 — P±p. Using the lemma of Ref (2j to remove the 
center-of-mass momentum term P, the reduced two-body Hamiltonian h becomes 
'Hn/2- If we now apply this simple bound to the gravitational problem, V{r) = —v/r, 
then we find, from Q and the one-body lower bound |2Jl, 



E > Nm 




(AT - 1)2^2 



{N - l)v < 1. 



(5) 



4. A^-particle problems in terms of Jacobi relative coordinates 

Our Hamiltonian H does not have the kinetic energy of the A^-particle system's center- 
of-mass removed. Thus, its eigenvectors are subject to two fundamental symmetries: 
translation invariance, and boson permutation symmetry (in the individual-particle 
coordinates {ri, r2, . . . , rjv}- Jacobi relative coordinates may be defined with the aid 
of an orthogonal matrix B transforming old ({r^}) to new {{Pi}) coordinates by 
[p] = B[r]. The first row of B, with all entries Bu — 1/\/N, defines a center-of- 
mass variable pi, its second a pair distance p2 = (ri — r2)/\/2, and the fcth row 
{k > 2) first has k—1 entries B^i = 1/ \/k{k — 1), the fcth entry Bkk — —\/(k— l)/fc, 
and zero for all remaining entries. The momenta {tTj} conjugate to the {p,} 
read [tt] = (i?^^)"'^[p] = -B[p]. Now, for an A^-boson problem with an attractive 
potential V{r), let 5'(p2, P3, ■ • • , Pn) be the (still to be found) normahzed ground-state 
eigenfunction corresponding to the lowest energy E. Boson symmetry is a powerful 
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constraint that greatly reduces the complexity of this problem. Although a non- 
Gaussian wave function is not necessarily symmetric in the Jacobi coordinates, we do 
have the remarkable A'^-representability expressions [5] (Appendix A here) for i,j > 1 

{^,{p,-p^)^) ^S.,{^,pI^) (6) 

(tT, • TTj)*) =S,j{^,7vl^). (7) 

We now choose to use pat and pn-i and we obtain the reduction 



E = {^,H<Sf) = ^ 



N 



N 



y \/Pw + "^^ + y VPat-i +m^+ lV{\rN^i - rjv|) 



where 7 — ^N{N — 1). We first note that the last row of B (which defines pn) is 

given by 

[a,a,a, . . . ,a,—{N — l)a], where a = l/y^N{N — 1). 

Now we express the equation for E in terms of new coordinates defined by the following 
relations 



R 

r 



P a 
a -P 



Pn 

Pn-1 



P 
P 



where 



N 



N - 1 



> 1, /? 



N -2 
N - 1 



1 
2 

< 1, 5 



Consequently 



and 



2, + 0^ 



{N - l)a = a 



P 



Na = a, S ^ —, 1 
a 



r = |r|, 




P a 




a -13 




N 


2 






N -1 




N ' 




1 + = 


2A 



TTJV-I 



< 1. 



We note that because of the boson symmetry, and after removal of tti, we have 
generally that 

(/(|PA'-i|)> = (/(|P - SP\)) = (/(|p + 6P\)) = (/(Ip^D), 

where f{p) is any appropriate kinetic-energy function. The expression for the energy 
now has the form 



E 



or, equivalently. 



6P) 



N 



{p- spy +"fV{r) ) , 



E^{n), where Ti = {p + 5Py + + -/Vir). (8) 

The Hamiltonian H is bounded below by the simple bound 0.^/2 of Eq. This 
result is proved in Appendix B. We now consider the eigenequation Hip — £ip, where 
'i/'(r,R) retains some of the boson symmetry implications of the full wave function. 
What we need are immediate consequences of Eq. namely (p^) — (P^), (p-P) ~ 0, 
and hence 



((p + <5P)2) = (p2 + s'p^) = (1 + s^){p^) = 2A(p2). 



(9) 



Gravitating semirelativistic N-boson systems 



6 



5. Improved lower energy bound 

The lower bound that improves the N/2 bound @ and forms our main resuh is 
specifically for the attractive pair potential V{r) — —v/r, v > 0. Suppose that ?/)(r, R) 
is the exact lowest eigenfunction of the Hamiltonian ?i in (|SJ). We look in a restricted 
domain "D in the Hilbert space L^{R^), one that keeps some of the original boson 
symmetry, namely, we assume Eq. This allows us to apply the same reasoning 
as we did for the one-body problem because the kinetic term is squared and thus the 
final form of the expectation values involves only the conjugate variables {r, p}. Let £ 
be the minimum of {ijjj'Hip) corresponding to normalized ijj satisfying £ is a lower 
bound to E, £ < E, because, if 4' is the exact normalized A^-body wave function, then 

E = i?*) ^ H*) > T^V) = 
The eigenvalue equation of 7i, Tiip — £'ip, explicitly reads 



Squaring these equal vectors and use of lO, and also the scaling change {r,p} 
{2r,p/2}, imphes 



By comparing this with the one-body case |(5J), we obtain 



E>£> Nm]J ^ ^ — , 71)2 < 1. ^^Q^ 

Figure Q illustrates the improvement by the lower bound (|10|) over the previously 
available simple lower bound Q. 

6. Gaussian upper energy bound 

In order to find an upper bound, we follow Ref. |1(J) and use a Gaussian wave function, 
which we write in the form 



where a > 0, while C guarantees the normalization of $. The boson symmetry of the 
trial function allows us to write E < Eg = H^) , where we have 



EG^[<i>, 



A^Jp2, + m2+7l/(|ri-r2|) 



$ (11) 



NJX7T% + m^+-fV{\V2p^\) 



$ . (12) 



The symmetry of the Gaussian function in the relative coordinates and the 
factoring property allow us to replace ttjv by 772 = p, and, setting r = we find 
explicitly 
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V 



Figure 1. Energy bounds E{v) (in dimcnsionless units) for tlie semirelativistic 
A'^-body problem with N = 5 and mass m = 1, as functions of the dimcnsionless 
coupling parameter v. The figure shows the simple lower bound (SL), the improved 
lower bound (L), and a scale-optimized Gaussian upper bound (U). 



where 0(r) = (a/Tr)'^/^ exp(— ^ar^). The kinetic-energy expectation value may be 
expressed in terms of the modified Bessel function of the second kind Ki{x) |26|. 
Evaluating the integral, for convenience, in momentum space, yields 

where fj, = m{2N/[{N — while g{x) is defined by 

/ + 00 
dtt^y/ix + t"^ exp(-^2). 
-oo 

The potential-energy expectation value reads for the case of a gravitational pair- 
interaction potential V{r) = —v/r: 



Nmv 12 
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With the expectation values in this form, 
variational parameter. We arrive at 



we may regard the parameter /l( as a 



E < Nm 



2 . 

mm 

TT fJ.>Q 



N >2. 



(13) 



The analytical fact that the function [Zglp^ /A) — a]/ fi has a minimum only if a < 2 
entails the constraint ti < 2/^. Figuren]confronts the two lower bounds jSj and ()10|l 
with this "scale-optimized Gaussian" variational upper bound. 




0.2 0.4 0.6 0.^ 



1.2 1.4 



Figure 2. Upper and lower bounds to the ground-state energy E{c) (in 
dimensionless units) of the gravitational A^-particle problem, N = 2, 3, . . . , 6, 
with mass rrt = 1, as functions of the dimensionless coupling parameter c. 



7. Results and conclusion 

Any consideration of an arbitrarily large number, TV, of self-gravitating bosons, that 
is, the inspection of the limit N oo, obliges us to force the gravitational coupling 
V to decrease for increasing N in some appropriate manner, in order to avoid the 
( "relativistic" ) collapse E — oo for N oo. Keeping v fixed, inevitably implies 
collapse [221 ■ Consequently, following Lieb and Yau ^3], let v diminish, with increasing 



Gravitating semirelativistic N-boson systems 



9 



according to v = c/N, where c is some constant gravitational interaction parameter. 
In this case the term {N — in the simple lower bound ||SJ) will be replaced by A^c^ 
(with c < 1), the improved lower bound (|10|l is modified by replacing 71;^ by ^Ac^ 
(with c < a/2) j and in the Gaussian upper bound H13I) we must substitute yj\/2c (with 
c < 2\/2) for ^J^v. Figure |21 shows resulting energy curves E{c) for N = 2, 3, . . . , 6 
and m = 1. For the lower bounds we require for the dimensionless parameter c < \/2; 
the upper bounds are valid, for all iV, if c < 2^/2. 

Some special cases are of interest. For small coupling c, both bounds are parabolic 
in shape. Explicitly, they read 

/ Ac^\ / Ac^ 

Nm 1 < E(c) < Nm 1 

\ 16 / \ 67r 

These energy bounds coincide with previous findings for the corresponding 
nonrelativistic problem; cf. Eq. (4.2) of Ref. ^S]- Meanwhile, for all couplings c < \/2 
(such that both constraints on c are satisfied), in the limit iV t 00 we have A — > 1, 
and the general results and 113|) provide bounds immediately to E{c)/Nm for 
arbitrarily large N. 

The energy spectrum of a system of nonrelativistic identical bosons experiencing 
any attractive pair interaction is bounded from below by the lowest energy of a 
specially scaled one-body problem [2]. For the harmonic oscillator, this bound yields 
the exact energy; such results are facilitated by use of (orthogonal) relative coordinates. 
Tight energy bounds have recently been obtained jlO| for semirelativistic problems 
with oscillator pair interactions. In this analysis we have constructed both lower and 
upper bounds to the lowest energy of a semirelativistic system of identical gravitating 
bosons. Our bounds reduce to their nonrelativistic counterparts in the limit of weak 
coupling. 
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Appendix A 

For definiteness, we consider the relative momenta tt^ and tt^ with i,j > 1. Since 
[tt] — B[p\, we have 

\ k / \ k / k k^l 

By using the boson symmetry of the wave function we have 



\ k J \k=^l J 



(A2) 
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But orthogonality of the rows of B to the first row tells us 

k \ k / \ I / \ k J \kM 



B,i I . (A3) 



If i = j, then the orthogonality of the matrix B tells us that 

Y,Bl = l = -Y,B,kBu. 

k k^l 

Meanwhile, if i ^ j, we know that the corresponding rows of B are orthogonal; hence, 
in this case, Bn^Bj^ — 0; thus, both of the coefficients in (A2) are zero. Hence we 
conclude for all i, j > 1 

{-Ki ■ TTj) = SijiTzl) ^ 6ij {(pI) - (pi • P2)) . (A4) 

The proof for the relative coordinates {p^} is identical. 

Appendix B 

The vectors P and p (with tti removed) define a plane; we let k be a unit vector 
perpendicular to this plane. Then, for example, we have 

(to^ + (p + ,5P)2)^ ^ |mk + p + (5P| (Bl) 

In this notation the expectation of the kinetic-energy term in 7Y, defined in Eq. 
reads 

(|mk + p + ^P|) = (|mk + p-(5P|) . (B2) 

Now, we may write 

mk + p = i[(mk + p + (5P) + (TOk + p-(5P)]. (B3) 
The triangle inequality then tells us 

Imk + p| < i |mk + p + 5P\ + ^ |mk + p - 5P\ . (B4) 

If we now look at mean values, we see from (B2) that 

(|mk + p|) < (|mk + p + 5P|), (B5) 

or, equivalently 

(^{m^+p^y)<(^{m^ + {p + SPfyy (B6) 

Thus we conclude that {H) = {H) > (7Y7v/2)i where the reduced one-body Hamiltonian 
Wiv/2 is given by 

nN/2^N{m^+p''y +jVir). (B7) 
We note the spectral inequality H > 7i7v/2 with equality only when N = 2. Thus 

E = (H) = (n)>{nN/2)- 
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